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ON  THE  CAUSE  OF  X-LINE  FORMATION  IN  THE  NEAR-EARTH  PLASMA  SHEET: 


RESULTS  OF  ADIABATIC  CONVECTION  OF  PLASMA-SHEET  PLASMA 


G.  M.  Erickson 


Department  of  Space  Physics  and  Astronomy 
Rice  University,  Houston,  TX  77251 

Abstract*  Self-consistent,  static-equilibrium  solutions  are  presented 
for  two-dimensional  magnetospherlc-magnetic-field  configurations  with  iso¬ 
tropic  thermal  pressure.  These  solutions  Include  a  dipole  field  and  are  not 
restricted  to  the  asymptotic  theory.  Adiabatic  convection  of  plasma  sheet 
flux  tubes  is  modeled  as  a  series  of  static-equilibrium  solutions  in  which 
flux  tubes  conserve  their  PVY  as  they  convect,  which  resulted  in  time  depen¬ 
dent  magnetospheric  configurations.  Specifically,  it  is  found  that  a  deep 
minimum  in  the  equatorial  Bz  develops  in  the  inner  plasma  sheet,  thereby 
causing  the  magnetic-field  configuration  to  become  more  stretched  and  tail¬ 
like  in  time.  These  results  suggest  X-line  formation  in  the  inner  plasma 
sheet  as  a  consequence  of  lossless,  adiabatic  convection  of  plasma  sheet  flux 
tubes. 
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Introduction 


X-line  formation  is  now  regarded  as  one  of  the  primary  candidates  for  the 
basic  physical  process  to  account  for  the  observed  energy  dissipation  associ¬ 
ated  with  magnetospheric  substorms.  It  is  clear  that  changes  in  the  solar  wind 
can  trigger  magnetospheric  substorms  (see,  e.g. ,  Akasofu,  1980].  However,  it 
is  not  clear  that  magnetospheric  substorms  are  always  associated  with  specific 
changes  in  the  solar  wind.  For  this  reason  a  substorm  mechanism  which  does 
not  necessarily  require  an  external  trigger  for  X-line  formation  is  attrac¬ 
tive.  Schindler  [1974]  offered  the  following  scenario  for  the  substorm  mecha¬ 
nism.  During  the  substorm  growth  phase,  free  energy  is  accumulated  in  the 
tail,  and  the  tail  becomes  more  and  more  unstable  to  perturbations  that  try  to 
create  neutral  lines.  The  presence  of  a  normal  magnetic  field  component  in 
the  current  sheet  inhibits  the  instability.  But  when  the  plasma-sheet  becomes 
sufficiently  thin  and/or  the  normal  magnetic  field  component  becomes  suffi¬ 
ciently  small,  some  breakup  mechanism,  such  as  the  ion-tearing  mode,  becomes 
operative  leading  to  X-line  formation  and  substorm  onset  [see  also  Nlshlda  and 
Nagayama,  1973;  Russell  and  McPherron,  1973;  Hones,  1977].  Erickson  and  Wolf 
[1980]  have  presented  the  argument  that  approximately  lossless,  adiabatic, 
earthward  convection  of  plasma-sheet  plasma  on  closed  field  lines  is  necessar¬ 
ily  time  dependent.  They  hypothesize  that  if  earthward  convection  occurs  in 
the  tail,  even  if  solar  wind  conditions  were  steady,  magnetotail  field  lines 
would  become  more  and  more  stretched  and  tail-like  in  time  resembling  the 
growth  phase  scenario  above.  Schindler  and  Birn  [1982]  and  Birn  and  Schindler 
[1983]  have  self-consistently  modeled  the  quasi-static  evolution  of  tail-like 
configurations  within  the  framework  of  the  asymptotic  theory.  They  found  that 
in  the  absence  of  unrealistic  plasma  loss,  earthward  convection  is  time  depen¬ 
dent  and  drives  the  tail  toward  instability. 

As  discussed  by  Erickson  and  Wolf  [1980]  it  is  the  presence  of  the 
earth's  dipole  field  which  prevents  approximately  lossless,  adiabatic  con¬ 
vection  of  plasma-sheet  flux  tubes  from  proceeding  in  a  time-independent 
manner.  Specifically,  it  is  the  fact  that  plasma  cannot  expand  very  far  along 
a  line  of  force  in  the  presence  of  a  dipole  field  which  causes  the  tall  con¬ 
figuration  to  become  more  stretched  and  tall-like  as  convection  proceeds.  The 
asymptotic  theory  excludes  some  (B#V)S  magnetic  tension  terms  in  the  force 
balance  of  the  system  and,  therefore,  cannot  accurately  represent  the  earth's 
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dipole  field  and  its  critical  role  as  the  endpoint  of  plasma-sheet  flux  tubes. 
In  this  paper,  preliminary  results  are  presented  of  self-consistent  modeling 
of  quasi-static  convection  of  plasma-sheet  flux  tubes  in  two  dimensions. 
These  models  include  a  dipole  field  and  are  not  restricted  to  the  asymptotic 
theory.  The  results  of  this  modeling  confirm  the  results  of  the  previous 
efforts  mentioned  above.  Also,  the  results  lend  support  to  the  Schindler 
[1974]  growth  phase  scenario,  although  they  do  not  exclude  the  possibility  of 
other  substorm  scenarios  [e.g.,  Coroniti  and  Kennel,  1972;  Akasofu,  1980; 
Atkinson,  1980]. 


Quasi-Static  Convection 


To  model  convection  of  plasma-sheet  flux  tubes  we  make  the  following 
assumptions:  (l)  plasma-sheet  ions  are  in  bounce  equilibrium,  (2)  thermal 

plasma  pressure  is  isotropic,  (3)  inertial  forces  are  small  compared  to  pres¬ 
sure  gradients,  and  (4)  convection  is  lossless.  Under  these  assumptions  we 
seek  static-equilibrium  solutions  of  the  momentum  equation  and  Maxwell's 
equations, 

3  x  g  -  VP  (1) 

J-f-Vxg,  (2) 

0 

and 


V  •  g  -  0.  (3) 

♦  a 

In  two  dimensions,  the  magnetic  vector  potential  is  A  -  Aey  (GSM  coordinates 
in  units  of  earth  radii  are  used  throughout)  with  Bx  ■  -3A/3z,  Bz  -  3A/3x,  and 
( 1 )— C 3)  are  rewritten  as 


V2A 


-v 


o 


dP£Al_ 

dA 


6(x)5(z), 


(4) 


with 


P, A  ■  constant  along  field  line 


In  (4),  the  delta  function  represents  the  dipole  source*  Numerical  solutions 
of  (4)  are  obtained  on  a  grid  displaced  from  the  origin,  and  A  is  decomposed 
into  its  dipole,  A<j,  and  plasma  current  source,  A j ,  parts.  Hence,  we  actually 
find  solutions  of 


»*A.  -  -u  ^ 


Adjustable  parameters  of  these  models  Include  the  equatorial  plasma  pressure, 
the  thickness  of  the  plasma-sheet,  and  the  location  of  the  magnetopause. 

By  the  frozen-in-flux  theorem  and  the  fact  that  in  two  dimensions  both  P 
and  A  are  constant  along  a  magnetic  field  line,  we  can  use  A  to  tag  the  plasma 
as  it  convects.  Also,  under  the  assumptions,  the  plasma  convects  so  as  to 
keep  its  PV^  constant,  where  y  is  the  adiabatic  index  (which  we  chose  to  be 
2),  and 

V  -  Jft  ds/B  (6) 

is  the  volume  of  a  flux  tube  of  unit  magnetic  flux.  Thus,  convection  is 
modeled  as  a  time-sequence  of  static-equilibrium  solutions  of  (5)  such  that 
each  solution  shares  the  same  PVY(A).  That  is,  we  require  that 

PV^(A)  ■  independent  of  time.  (7) 


Results 

Some  examples  of  self-consistent,  two-dimensional  ma gnet os pheric -magnet¬ 
ic- fie  Id  configurations  are  shown  in  Figure  1  for  different  pressure  distribu¬ 
tions,  plasma-sheet  heights,  and  magnetopause  locations.  Figures  2-5  show  the 
equatorial  pressure  distribution,  the  equatorial  magnetic  field,  PV2(x),  and 
PV2(A),  respectively,  for  these  models.  Note  that  in  order  for  approximately 
lossless,  adiabatic  convection  to  occur  in  some  region  of  the  tail  in  a  time- 
independent  manner,  PVY (x)  would  have  to  be  constant  in  that  region.  However, 
for  static-equilibrium  configurations  that  satisfy  (1)  and  include  a  dipole 
field,  PV^(x)  (for  y  >  5/3)  is  not  constant  but  increases  with  distance  down 
tail.  Thus,  as  earthward  convection  proceeds,  flux  tubes  entering  some  region 


pro 


would  have  a  larger  particle  content  than  did  their  predecessors  in  the 
region.  Conditions  (1)  and  (7)  would  then  require  that  the  magnet ospheric 
configuration  vary  in  time. 

Starting  with  some  self-consistent  solution  of  (1)  or  (5)  at  t  *  0,  an 
electric  field  Emp(t)  is  imposed  at  the  magnetopause  allowing  magnetic  flux  to 
enter  the  magnetopause  and  forcing  earthward  convection  of  plasma-sheet  flux 
tubes.  In  practice,  the  value  of  A  at  the  magnetopause,  A,^,  is  changed,  the 
equatorial  plasma  pressure  is  adjusted,  and  new  solutions  of  (5)  are  obtained 
until  (5)  and  (7)  are  satisfied  simultaneously.  Thus  a  time  sequence  of 
static-equilibrium  solutions  of  (l)-(3)  is  constructed  with  time  parameterized 
by  the  value  of  Amp.  Forcing  Models  A  and  B  to  convect  in  this  manner  result¬ 
ed  in  Models  A'  and  B'  shown  in  Figures  6  and  9.  In  Model  B'  modeling  was 
performed  only  out  to  x  *  -60.5  Rg  where  the  height  of  the  plasma-sheet  was 
chosen  so  that  PV2(x)  had  zero  slope  there  at  t  *  0.  Beyond  x  *  -60.5  Rg,  PV2 
was  assumed  constant.  For  t  >  0,  the  height  of  the  plasma-sheet  was  chosen 
such  that  the  amount  of  magnetic  flux  contained  between  the  top  of  the  plasma- 
sheet  and  the  magnetopause  was  constant  in  time.  For  quiet  times  and  a  con¬ 
stant  electric  field  of  10“4  V/m  at  the  magnetopause,  the  three  "snapshots”  of 
Figures  6  and  9  for  Amp  -  0,  -12,  and  -22  correspond  to  t  *  0,  1.5,  and  2.7 
hrs,  respectively.  For  an  electric  field  of  5  x  10“4  V/m  the  respective  times 
would  be  t  ■  0,  0.3,  and  0.54  hrs.  The  conversion  of  the  other  physical  quan¬ 
tities  from  the  units  used  to  SI  units  is  given  in  the  Appendix. 

As  expected,  the  adiabatic  convection  of  plasma-sheet  flux  tubes  resulted 
in  time-dependent  configurations,  consistent  with  the  earlier  conclusions  of 
Erickson  and  Wolf  [1980]  and  Schindler  and  Bim  [1982]  that  steady-state  con¬ 
vection  is  theoretically  unlikely  in  the  magnetotail.  Being  unable  to  convect 
in  a  steady  state,  the  field  strength  in  the  lobe  increases  as  flux  is  piled 
up,  and  the  configuration  becomes  more  tail-like  in  time.  Figures  7  and  10 
show  the  evolution  of  the  equatorial  pressure  distribution  for  Models  A'  and 
B'.  The  drift  of  plasma  around  the  earth  was  accounted  for  by  defining  the 
equatorial  pressure  at  the  origin  to  be  the  maximum  pressure  of  the  system  for 
all  times.  Thus  the  pressure  was  allowed  to  reach  this  value  but  not  exceed 
it.  Figures  8  and  11  show  the  evolution  of  the  equatorial  magnetic  field  Be 
for  Models  A'  and  B*.  Note  the  minimum  in  Be  that  develops  (and  gets  deeper 
with  time)  near  -10  Rg,  corresponding  to  the  stretching  of  the  field  in  the 


near-earth  part  of  the  plasma-sheet  as  higher  and  higher  content  flux  tubes 
are  convected  into  the  region.  Also,  the  plasma  current  density  has  approxi¬ 
mately  doubled  In  this  region  from  Amp  »  0  to  A^  -  -22.  It  Is  In  this  near- 
earth  plasma-sheet  region  that  the  tall  appears  least  stable. 

Admittedly,  the  manner  in  which  the  drift  of  plasma  around  the  earth  Is 
treated  does  not  seem  very  realistic.  This  feature  was  chosen  merely  for 
numerical  simplicity  and  convenience.  Condition  (7)  was  strictly  enforced 
tailward  of  the  P  =  constant  region,  while  it  was  not  enforced  inside  the  P  = 
constant  region.  This  resulted  in  a  sharp  inner  edge  of  the  current  sheet 
with  current  density  increasing  from  zero  to  its  peak  value  within  only  a 
few  Rg.  The  westward  currents  in  the  plasma-sheet  contribute  a  tk  tive 
Be  earthward  of  the  peak  in  the  current  density.  As  higher  and  highei  ontent 
flux  tubes  convected  into  the  near-earth  part  of  the  plasma-sheet,  west¬ 
ward  current  increased  resulting  in  a  local  Be  minimum  just  earthwa  the 
current  density  maximum.  In  a  more  realistic  treatment  the  drift  ot  plasma 
out  of  the  noon-midnight  meridian  plane  and  around  the  earth  would  still  occur 
within  only  a  few  Rg,  but  the  artificial  constraint  that  the  pressure  not 
exceed  a  certain  constant  value  would  be  removed.  In  that  case  a  peak  in  the 
plasma  pressure  might  develop  just  tailward  of  the  inner  edge  of  the  plasma- 
sheet.  The  current,  -dP/dA,  would  change  sign  at  this  peak  in  a  closed-field- 
line  configuration.  As  higher  and  higher  content  flux  tubes  convect  into  the 
region,  the  westward  current  tailward  of  the  peak  and  the  eastward  current 
within  the  inner  edge  would  increase,  resulting  in  a  Be  minimum  near  the 
plasma  pressure  peak.  Thus,  we  would  expect  that  a  more  realistic  treatment 
would  not  qualitatively  change  the  results.  However,  details  such  as  the 
exact  location  or  width  of  the  Be  minimum  might  be  different. 

A  rectangular  magnetopause  was  also  chosen  for  numerical  simplicity. 
Results  of  the  static  modeling  show  that  the  choice  of  dayslde  magnetopause 
has  little  effect  on  the  tail  configuration.  Also,  flaring  of  the  tail  magne¬ 
topause  affects  the  PV^(A)  curves  in  much  the  same  way  as  the  height  of  the 
plasma-sheet  does.  Thus  we  expect  that  the  general  behavior  of  the  convection 
models  presented  here  will  also  occur  for  more  realistic  magnetopause  shapes. 


Summary 


The  effect  of  slow,  sunward  convection  In  a  magnetospherlc  plasma-sheet 
has  been  Investigated  using  computed  two-dimensional,  force-balanced  magnetic 
field  configurations,  Including  the  earth’s  dipole  field.  The  results  confirm 
the  earlier  conclusions  of  Erickson  and  Wolf  [1980]  and  Schindler  and  Blm 
[1982]  that  approximately  lossless,  adiabatic  convection  of  plasma-sheet  flux 
tubes  is  a  time-dependent  process.  In  this  process  the  magnetotall  becomes 
more  stretched  and  tail-like  as  convection  proceeds,  resulting  In  configura¬ 
tions  more  unstable  to  perturbations  which  try  to  create  neutral  lines.  The 
results  support  the  Schindler  [1974]  growth— phase  scenario,  although  they  do 
not  exclude  other  possibilities.  Preliminary  results  suggest  that  forcing 
plasma— sheet  plasma  to  adiabatlcally  convect  sunward  leads  to  a  buildup  of 
magnetic  energy  in  the  magnetotall.  A  minimum  In  the  equatorial  field 
strength  develops  in  the  near-earth  part  of  the  plasma-sheet,  and  that  minimum 
rapidly  deepens  as  convection  proceeds,  which  obviously  suggests  formation  of 
a  near-earth  X-line. 


Appendix 

The  unit  of  distance  used  here  is  Rg  (6.38  x  106  m).  The  other  physical 
quantities  are  given  in  arbitrary  units  which  can  be  converted  to  MKS  units  as 
follows: 


P  -  6P/(y0RE*t), 

A  -*■  d^A/Rg, 

B  d^B/Rg2, 

V  V/(dl/2Rg3). 

Time  is  parameterized  as  the  value  of  A,^,  the  value  of  A  on  the  magnetopause. 
Given  the  electric  field  at  the  magnetopause,  Eap(t),  time  (in  seconds)  is 
determined  from 
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Figure  Captions 


Some  examples  of  self-consistent,  two-dimensional,  magnet ospheric- 
magnetic-f ield  configurations.  Tail  field  lines  are  shown  containing 
equal  amounts  of  magnetic  flux.  The  amount  of  magntic  flux  passing 
through  the  right  boundary  is  the  same  in  each  example.  (a)  Model 
A.  Using  the  appropriate  boundary  conditions  (P  and  A  at  the  magne¬ 
topause,  P  in  the  equatorial  plane,  Bx  (z  ■  0)  ■  0,  and  one-dimen¬ 
sional  force  balance  at  the  right  boundary)  a  modified  (P  *  constant 
for  0  >  x  >  -4.5  Rg)  Fuchs-Voigt  [1979]  model  with  k  «*  1.54  is  ob¬ 
tained.  The  Fuchs-Voigt  anlaytic  models  feature  an  exponential 
decline  of  the  physical  quantities  down  the  tail,  (b)  Model  B.  This 
model  is  the  same  as  Model  A  except  that  the  height  of  the  plasma- 
sheet  (dotted  line)  is  chosen  as  6  Rg  at  x  “  -60.5  Rg.  (c)  Model 
C.  In  this  model  the  equatorial  pressure  declines  as  |x|”1*2,  based 
on  the  observations  of  Behannon  [1968].  (d)  Model  D.  This  model  has 

the  same  equatorial  pressure  as  Model  A,  but  the  tall  magnetopause  is 
flared  (5.7°),  and  the  dayside  magnetopause  is  rounded. 

The  equatorial  pressure  profiles  for  Models  A-D.  Note  that  dP/dx  for 
Models  A,  B,  and  D  is  the  same. 

The  equatorial  magnetic  field,  Be,  for  Models  A-D. 

PV2  vs.  x  for  Models  A-D.  A  constant  PV^(x)  is  required  for  steady, 
lossless,  adiabatic  convection. 

PV2  vs.  A  for  Models  A-D.  The  approach  of  the  top  of  the  plasma- 
sheet  to  the  equtorlal  plane  implies  a  neutral  line  tallward  of  the 
right  boundary  in  Model  B.  When  Model  B  is  convected,  it  is  assumed 
that  flux  tubes  of  constant  PV2  enter  the  right  boundary  represented 
by  the  dashed-line  continuation  of  the  PV2(A)  curve. 
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Fig.  6.  Model  A'  showing  che  results  of  convecting  Model  A  under  the  con¬ 
straint  that  PV2(A)  remains  unchanged  (as  explained  in  the  text). 
The  two-dimensional  magnetic  field  configuration  if  shown  for  Amp  «  0 
(top  panel),  Amp  *  -12  (middle  panel),  and  Amp  -  -22  (bottom  panel). 
Tail  field  lines  are  shown  containing  equal  amounts  of  magnetic  flux. 
The  A  ■  0  field  line  is  labeled,  and  the  dashed  field  line  is  A  ■  30 
for  reference. 

Fig.  7.  The  equatorial  plasma  pressure  for  Model  A'  for  Amp  *  0,  -12,  and 
-22.  The  end  of  the  P  -  constant  region  is  at  x  •  -8.5  Rg  for 

A|qP  <  —12. 

Fig.  8.  The  equatorial  magnetic  field  for  Model  A'  for  Amp  -  0,  -12,  and  -22. 

Fig.  9.  Model  B'  shows  the  results  of  convecting  Model  B  in  the  same  manner 
as  Model  A  with  the  additional  condition  that  the  total  flux  contain¬ 
ed  between  the  top  of  the  plasma-sheet  (dotted  line)  and  the  magneto¬ 
pause  is  constant.  The  snapshots  shown  are  displayed  in  the  same 
manner  as  Figure  6. 

Fig.  10.  Same  as  Figure  7  for  Model  B'. 


Fig.  11.  Same  as  Figure  8  for  Model  B* 
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